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Abstract 

Title of Thesis: Heat Transfer Characteristics of Hypersonic Waveriders 

with an Emphasis on the Leading Edge Effects 

Name of Degree Candidate: Denis O. Vanmol 

Degree and Year: Master of Science, 1991 

Thesis Directed by: Dr. John D. Anderson, Jr. 

Professor 

Department of Aerospace Engineering 

The heat transfer characteristics in surface radiative equilibrium and the aerody- 
namic performance of blunted hypersonic waveriders are studied along two constant 
dynamic pressure trajectories for four different Mach numbers. 

In terms of performance, the inviscid leading edge drag was found to be a small 
(4 to 8%) but non-negligible fraction of the inviscid drag of the vehicle. Although the 
viscous drag at the leading edge can be neglected, the presence of the leading edge 
will influence the transition pattern of the upper and lower surfaces and therefore 
affect the viscous drag of the entire vehicle. 

Finally, for an application similar to the NASP, the present study demonstrates 
that the waverider remains a valuable concept at high Mach number if a state-of- 
the-art active cooling device is used along the leading edge. At low Mach number 
(<5), the study shows that surface radiative cooling might be sufficient. In all cases, 
radiative cooling is sufficient for the upper and lower surfaces of the vehicle if ceramic 
composites is used as thermal protection. 
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Chapter 1 INTRODUCTION 


In 1989, a panel of experts in hypersonic technology published a report 1 defining 
the most critical technologies of hypersonic flight for an application such as the 
NASP. They judged these to be : 

• the supersonic combustion ramjet engine 

• the technology of integrating such engine with the airframe 

• the structural concepts and high temperature low density materials required to 
achieve the desired weight fractions 

• guidance and control of such large, flexible, complex vehicle 

The present research, focused on the heat transfer evaluation of an hypersonic 
waverider aiming at the same mission as the NASP, is directly related to the third 
technological key through the definition of the thermal protection of the vehicle. 

Previous studies at University of Maryland^’ 4 have already demonstrated that 
viscous-optimized waveriders can be designed with lift-to-drag ratios higher than 
other conventional hypersonic configurations. Moreover, the waverider generates a 
relatively uniform flow on the lower surface making this concept very suitable from 
the engine integration point of view"’ . In another study** conducted at University of 
Maryland, results indicate that the point-design waveriders do not suffer any marked 
off-design problems; they indeed behave much like delta wings in such a cases. 

For all these reasons, the waverider concept is moving from an academic concept 
to a very practical one for consideration as hypersonic cruiser, hypersonic accelerator 
or aero-assisted vehicle for space mission. 

The present work is part of the continuing study of waveriders in the Hypersonic 
group at the University of Maryland. This study takes another logical step and 
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examines the problem of aerodynamic heating for this class of vehicles with an 
emphasis on the leading edge. 

Section 1 Purpose of study 

The purpose of this study is to demonstrate that the waverider concept remains 
a valuable option when the problem of aerodynamic heating is taken into account. 
This problem is probably more critical at the leading edge of a waverider than for 
any other vehicle. 

The leading edge of the waverider being in theory infinitely sharp, it will lead to 
an infinite value of the temperature. There is thus an absolute need for blunting the 
leading edge in order to decrease the temperature. But, due to the small waverider 
thickness near the leading edge, the accommodation of a finite radius also means that 
the planform shape will have to be modified. In this case, if the planform shape is 
too strongly affected, the flow on the lower surface will loose its original nature and 
the performance of the waverider will be degraded. The name of the game is thus to 
implement, according to a given waverider geometry, a radius distribution all along 
the leading edge which is such that we will not exceed the maximum temperature 
allowed for a given material. 

When dealing with a finite value of the radius along the leading edge, another 
major phenomena has to be taken into account : the transition along the attachment 
line of the leading edge. This phenomenon, is relevant not only along the attachment 
line in the spanwise direction , but it also has an impact on the transition pattern in 
the streamwise direction. 

Finally, we will demonstrate that the waverider remains a valuable concept if a 
state of the art active cooling device is used at the leading edge. 

In order to achieve this work, a computer code originally written by McLaughlin 4 
was extensively modified. The heat transfer characteristics for surface radiative 
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equilibrium, as well as the aerodynamic performance of waveriders with blunted 
leading edges, arc studied along two constant dynamic pressure trajectories for four 
different Mach numbers. 
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Chapter 2 WAVERIDER CONSTRUCTION 


The computer code originally written by McLaughlin was modified to implement 
the leading edge radius distribution and the heat transfer evaluation on the vehicle. 
Since this study does not require the modification of the procedure used to generate 
viscous optimized waveriders, we will extensively refer, for this part, to the original 
work of McLaughlin 4 . 

To generate a waverider, two basic functions must be performed. First, the 
flowfield around a generating body (cone here) is calculated. Second, when the 
generating flowfield is known, the stream function must be evaluated in order to 
define the flow streamlines. Once these are known, it is possible to carve any number 
of different waveriders from a single generating flowfield. 

Section 1 Generating flowfield solutions 

Calorically perfect case 

The conical flowfield solution for the calorically perfect gas is generated using 
a standard Taylor-Maccoll formulation 7 . An ordinary differential equation is solved 
numerically for a specific cone angle using an inverse approach where the shock angle 
is guessed and the cone angle that supports the shock is calculated. 

For a given shock angle, we can calculate the conditions immediately behind the 
shock. These conditions are used as initial values and the Taylor-Maccoll equation is 
integrated using a fourth order Runge-Kutta method from the shock to the body. Once 
the proper solution is achieved, all the flow properties such as pressure, temperature, 
velocity can be defined. 
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Chemical equilibrium case 

The same inverse approach is used to generate the chemically reacting conical 
flowfield 8 . In this case a set of three coupled differential equations combined with 
equilibrium high-temperature properties of air 9 is integrated using the same fourth 
order Runge-Kutta method. This system of equations is a chemical equilibrium 
analogous to the calorically perfect Taylor-Maccoll equation for conical flow at zero 
angle of attack. 


Section 2 Generation of waverider geometry 

The waverider geometry can be entirely defined when its planform shape, coin- 
cident with the shock, is known. From each point defining the planform shape , a 
streamline can be traced through the conical generating flowfield and the resulting 
stream surface becomes the lower surface of the waverider. From the same points, 
freestream streamlines can be created and will determine the upper surface (fig.l). 

Note that since the upper surface of the waverider is a freestream stream surface, 
the flowfield depicted above this surface, and the generating cone, do not exist when 
the waverider is actually in flight. The entire flowfield above the upper surface is 
strictly used in the design of the vehicle, and only the flowfield below the waveriders 
lower surface actually exists at flight conditions. 

Due to the waverider symmetry, only half of the planform shape needs to be 
defined. Each planform shape is defined by the projection, on the generating conical 
shock, of a curve passing through five points in the x-y plane (fig.2). 

By symmetry, one point is forced to lie on the y-axis. Another one is forced to 
lie on the conical shock. These five points represent in fact eight degrees of freedom 
(10 coordinates — 2 constraints) which will be used by the optimization procedure 
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to determine the best geometry associated with a given figure of merit ( maximum 
L/D, maximum volume, ...). 

Section 3 Definition of the leading edge geometry 

Up to here, the generated waverider has still an infinitely sharp leading edge. 
Since the waverider geometry of the upper and lower surfaces are now known we can 
define a process to accommodate some bluntness to the leading edge. Let’s consider 
a side cross section of our waverider (fig.3) at a given spanwise (Y=constant) station. 
The flow is coming from left to right. Moving in the Z direction from the original 
leading edge point to the trailing edge, we can define a radius R at any Z station in 
function of the local waverider thickness and the local relative angle of attack ( 8 ) of 
the lower surface. The waverider length, in a Y=constant plane is defined as L(Y) 
and represents the distance from the original leading edge point to the trailing edge. 
In the symmetry plane, L(Y) equals the imposed waverider length and L(Y)=0 at the 
wing tip. Thus, L(Y) varies in between these two limits according to the vehicle 
planform shape. 

The planform shape of a half waverider being defined by 51 points, we are now 
able to associate to each of these points a finite radius if we specify a certain amount 
of carving. The carving distance is measured from the original leading edge point 
to the point where the waverider thickness is evaluated. This distance, C*L(Y), can 
be expressed in terms of a carving coefficient C held constant or varying with L(Y). 
The most satisfying results were obtained with a carving coefficient varying with the 
ratio L(Y)/L. The function used is of the type : 

= -™ th < = (I> 
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where a, b, m are positive coefficients adjusted to obtain a reasonable variation 
of C((). For all the cases presented later, m was set to 0.9 and the value of a, b 
were calculated in such a way that C(0)=50% and C(l)=2%. This increasing carving 
coefficient, from the centerline to the wing tip, combined with decreasing values 
of L(Y) produces a slowly decreasing carving distance which preserves relatively 
well the planform shape (fig.4). A typical radius distribution is presented in figure 
5 and shows a slowly decaying radius distribution except near the wing tip where 
it decreases rapidly. Since L(Y) and the waverider thickness are equal to zero at 
the wing tip, the radius at that particular point is arbitrarily set to the value of the 
previous one. Since the radius distribution is not guaranteed to be a monotonically 
decreasing function we will enforce it. This is done by simply imposing that the 
value of R(Y) is always greater or equal than the value of R(Y+AY). Therefore, 
the carving distance previously defined can be seen for all cases as the maximum 
modification performed on the planform shape. 

Once the radius distribution is defined, we can also calculate the sweep angle 
distribution along the planform shape (fig-6). Note that when the sweep angle 
distribution is calculated on the original planform shape, the differences between 
the two distributions are marginal. The sweep angle at the nose is set to zero and the 
sweep angle of the last point (wing tip) is set to the value of the previous point. 

Section 4 Definition of the nose geometry 

As we have defined our leading edge up to now, the nose would be represented 
by the intersection of two swept cylinders. This configuration, obviously not realistic 
from a heat transfer point of view, was modified by assuming that an arc of circle 
(radius R x , center C x ) joined the stagnation point and the second point of the planform. 
This assumption leads us to consider the nose region as a 3D stagnation point defined 

by the two radii R x and R y (fig-7). 
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Therefore, we will always assume that the nose region will be situated between 
the two first points of the planform shape. The leading edge region starts at the 
second point and extends up to the wing tip. 
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Chapter 3 AERODYNAMIC ANALYSIS 
OF A WAVERIDER 


The aerodynamic analysis of a wavcrider requires both an inviscid and viscous 
analysis on the vehicle. 

For the inviscid part, the analysis is relatively straightforward since it only 
requires the integration of the pressure forces. The inviscid aerodynamic analysis 
of a waverider configuration involves the calculation of both the lift and drag forces 
due to the pressure on the surfaces of the body. Since this part has already been 
extensively explained 4 for the upper and lower surface, we will only consider here 
the calculation of lift and drag due to the leading edge. 

For the viscous part, the analysis is more complicated since it is necessary to 

obtain : 

• the transition front delimiting the laminar and the turbulent regions 

. the radiative equilibrium wall temperature associated to each point of the wa- 
verider surface 

• the shear stresses calculated at the corresponding wall temperature 

• the integration of the shear stresses on the complete vehicle 

Section 1 Leading edge inviscid lift and drag 

Half of the waverider planform is formed by 51 points where we already have 
defined a radius and a sweep angle (fig.8). But, to calculate the leading edge drag, 
the pressure distribution needs to be determined too. Since we are dealing with a 
hypersonic vehicle and a stagnation region, we know that the modified Newtonian 
theory does a good job of predicting the pressure distribution around blunt bodies. 
Moreover, it was also demonstrated that this pressure distribution is not affected 
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by the sweep angle and the flow nature. We can therefore state that the pressure 
distribution around the leading edge at a point i will vary as: 


Cp(r) = Cp a ( A,) cos J (r ) 


(2) 


where according to the modified Newtonian theory: 


C'Pa(Aj) = Cpt cos'(A,) 


(3) 


with 


C Pt = 


P02 ~ Poo 

Qoo 


TTierefore, at any station i, we will have : 


(4) 


Cp( T ) = Cpt cos ! (A ,) cos’(r) (5) 

The axial coefficient c" can be calculated from the integration of the pressure distri- 
bution in the reference axis x* , z as shown in figure 8 and 9. 

TE 

c; = i j (Cp,-Cp,)dz (6) 



Since 

Cp(r) - Cpt cos ? ( A , ) cos ; (t) (9) 
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we can replace it in the above equation and integrate over the boundaries. 


c 


n 

a. 


1 

3 


Cpt COS 2 A, 
[l + cos 6,) 


(2 + cos 6, (sin 2 6, + 2)) R, 


( 10 ) 


By definition, 


c 


ff 

a. 



floo^ 


qooRtH + costf.) 


( 11 ) 


Therefore, 

= ]-Cpt cos 2 A, (2 + cos 6, (sin 2 6, + 2)) (12) 

<7oo 3 

Since a" represents the axial force per unit span perpendicular to the planform shape 
at point i, only the component parallel to the x axis will participate to the drag of 
the vehicle. Defining 


floe 


— cos A, 

floo 


(13) 


where d" is the drag per unit span in the x direction. Assuming now that this quantity 
will vary linearly from d" to d'{ +] between s t and s, + , we can calculate the drag due 
to this segment of the leading edge. 


floo 


2</oo 


(14) 


Summing on all the segments, we obtain the total drag 


D_ 

floo 


E 


A 

floo 


(15) 


Finally, the leading edge drag coefficient for both sides is obtained as 

D 1 

GDlB — ^ n 

flOO Planform 


(16) 


it 



(17) 


n _ 
(‘Die ~ r ? 


Cpt 

3 s n 


an form 


[( /?, cos 3 A, (2 + cos 8 t (sin 2 8, + 2)) + 


/? 1+1 cos 3 A l+1 (2 + cos <5 I+ , (sin 2 £, +1 + 2)))(s l+1 - s,)] 

Due to the non symmetry of the leading edge design, this configuration will produce a 

small amount of negative lift. Applying the same assumption as above, this quantity 

can be expressed as : 


Cl,, e = 


C Pt 1 


3 5 


Planform 


imoi-1 

[(/2» cos 2 A, sin 3 8 t + 

U1 


(18) 


cos 2 A, +1 sin 3 6 1+1 )(s, + l - 5,)] 

In this simplified analysis, the calculation of the lift and the drag due to the leading 
edge has neglected the presence of the anhedral angle which is usually small on 
waveriders. This angle, when combined with the axial force a" would eventually 
produce a small amount of positive lift. 


Section 2 The prediction of transition 

Transition at the attachment line of the leading edge 

Let’s consider a simple cylinder of very large span placed in a uniform flow such 
that the flow is normal to the spanwise generators. In this situation, the boundaty 
layer which forms on the body will be two-dimensional and begin at the stagnation 
line. In this case, we would always expect to find a region of laminar flow near 
the stagnation line because the Reynolds number based upon local velocity tends 
to zero as the stagnation line is approach. If the cylinder is now slightly swept at 
an angle A (fig. 10), the flow outside the boundary layer becomes three dimensional. 
Therefore, the application of sweep produces an additional velocity component tangent 
to the cylinder and we must determine what effect sweep will have on the transition 
behavior of this flow. 

If we now assume that the Reynolds number is infinitely large, the boundary 
layer on the cylinder surface has zero thickness. In this case we can consider that a 
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particular stream surface will correspond to the swept cylinder surface. This stream 
surface will encounter first the line AA and then divides with one branch following 
the upper surface and the other following the lower branch of the body. The line AA 
is called the “attachment line”. In the 2D case, the attachment line AA is also the 
locus of two-dimensional stagnation points. In the case of a finite Reynolds number, 
the attachment line is taken to be the projection of the dividing stream surface onto the 
surface of the body (fig.l 1). In order to characterize the transition along the leading 
edge. Poll 11 defines H as a characteristic Reynolds number which can be used as a 
similarity parameter for the attachment line boundary layer even when the flow is 
transitional or turbulent. R is defined as follow : 

R=y^ (i9) 

v c 


and 


V = 


{du e /dx) gm0 J 


1/2 


( 20 ) 


is a boundary layer characteristic length. 

A large bulk of Wind-tunnel data and flight test experiments demonstrated that 
in the presence of large boundary layer tripping devices such as streamwise end- 


plates, boundary layer fences, isolated three-dimensional roughness elements or two- 
dimensional trip wires, transition at the attachment line begins when R exceeds 
approximately 250. 

When the surface on which the attachment line boundary layers form is smooth 
and there are no sources of upstream disturbances, the transition will be the result 
of the amplification of small fluctuations always present in the freestream. In such a 
case, the laminar flow will selectively amplify certain disturbances characterized by 
their frequency and wave number. 


Therefore, wave packets will appear in the boundary. As the wave packets are 
convected along the leading edge, the disturbance amplitudes increase and eventually 
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break down to form turbulent spots. The occurrence of turbulent spots, or bursts 
marks the onset of transition. This process is similar to that which is observed in 
two-dimensional flows with low freestream disturbances. For small tripping devices 
or smooth surfaces, it was found that transition was delayed until much higher values 
of R were reached ( R > 570 ). 

In order to examine more closely the effect of surface roughness or tripping 
devices. Poll extended the swept cylinder experiments to investigate the response of 
the laminar attachment line boundary layer to the presence of trip wires of various 
diameters. The results of this exercise are summarized in figure 12 which shows 
the value of R necessary for the observation of the onset of transition at a certain 
distance s from a trip wire with a diameter d. 

These results indicate that the transition behavior of the infinite swept attachment 
line may be separated into four distinct regimes : 

First regime : d A/ < 0.8 

For very small values of d, the results are indistinguishable from the free transition 
results. The values of R for transition onset are high and the turbulence is the result of 
an instability of the laminar flow to the small disturbances present in the freestream. 
This region is bounded by a critical value of d A/ and for values of d /77 greater than this 
limiting level, the transition location is determined by disturbances introduced by the 
wire. The behavior is qualitatively similar to that observed in the flat plate boundary 
layer where the concept of a maximum tolerable roughness height is well established. 
However, roughness height comparisons have indicated that the attachment line is 
much less sensitive to a trip wire. 

Second regime : 0.8 < d /77 - 1-6 

When the trip wire is sized such that transition occurs between an R of 400 
and 600, it is found that bursts of turbulence still occur first at large spanwise 
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distances downstream of the trip. However, the laminar perturbation which precedes 
the turbulent burst is no longer a packet of Tollmien-Schlichting waves but discreet 
disturbances traveling along the attachment line. Disturbances emanating from the 
trip wire dominate the transition process and the value of R for transition exhibits 
a strong dependence upon d/7/. 

Third regime : 1.6 < d/r/ < 2.0 

If the trip diameter is increased still further, the turbulent bursts originate at the 
trip wire and then propagate indefinitely along the span. The onset of transition is 
observed at all spanwise stations at the same value of R. 

Fourth regime : d/r; > 2.0 

In this case, the large trip wires or any gross disturbance like a wing-fuselage 
junction produce a similar response in the boundary layer. If the boundary layer is 
subjected to gross disturbances from an isolated source, transition begins simultane- 
ously at all those attachment line stations which are more than 20007/ downstream 
from the disturbance source as R exceeds 250. 

When the trip size is such that transition begins at an R of less than 250, the 
turbulent bursts still originate near the trip wire but decay and ultimately disappear as 
they convect along the attachment line if R stays below 250 at all those attachment 
line stations. Therefore, an R of 250 represents a lower limit for the indefinite 
propagation of turbulence along a swept attachment line when the boundary layer is 
subjected to a very large disturbance. This regime is the one of most direct practical 
interest for the following reasons: 

• when wing/fuselage or tail/fuselage or nacelle/fuselage junctions will be present 

on the vehicle, they will act as a gross disturbance and we will have to consider 

that the transition along the attachment line will take place in this regime. 
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if such a wing/fuselage junction is not present as in the case of the waverider, 
this regime will also be representative of the roughness state (d) of the leading 
edge on a conventional vehicle as we will see it later on. 


The variation of R with d/r; and s /77 for complete turbulence is extracted from reference 
12 and given in figure 13. By using the Emmons probability theory for transition, 
a model of the attachment line transition process was developed by Poll 13 . He 
concluded that, in the case of gross contamination, the intermittency at an infinite 
swept attachment line depends upon the parameters R, r\ and s in such a way that : 


P = I - EXP 


R 


- 24 5 \ / s 
233 / \Rrj ^ 


with -=— > 8 and R > 245 (21) 

Rrj 


Therefore, while under certain circumstances the onset of transition may be 
independent of spanwise position s, the completion of transition is always dependent 
upon s. 

Since the waverider presents no discontinuity in its planform shape, its attachment 
line along the leading edge is not subject to boundary layer contamination due to a 
wing fuselage junction. It was thus decided to consider in the frame of this study 
two test cases : 

• a so called “smooth” leading edge where d is such that d /77 is always less than 
0.4 (first regime). 

• a so called “rough” leading edge where d is such that d /77 is always higher than 
2.4 (fourth regime). 

Since the transition for the smooth leading edge cannot be described by the intermit- 
tency function F, it was decided for both cases to use the values of R and s /77 as 
given in figure 12 and 13 for the onset and the completion of turbulence. 
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The effect of compressibility 

Most of the results described up to here were conducted by Poll at subsonic speed. 
He later extended his work 13 to include the effect of compressibility on the transition 
at the attachment line. In his study, he only considered test cases in which transition 
was clearly the result of gross upstream contamination (fourth regime). He found 
that all the available data could be correlated by a single value of the attachment 
line similarity parameter R , provided that the temperature dependent properties were 
evaluated at some reference temperature T* (fig. 14). The compressible version of R 
was defined as R„ in the following way : 

«• = «/§) <22) 


where 


T, = T e + 0. 1 0{T W -T e ) + 0.60(7; - T e ) (23) 

* 

Later on, experiments conducted at Mach = 3.5 in the Pilot Low-Disturbance 
Wind-tunnel at NASA Langley 14 demonstrated that Poll’s criteria was valid not only 
in the fourth regime but also for the three others when compressibility effects were 
taken into account. 

The results show that end plates or large trips near the upstream end of the 
cylinder cause turbulent flow along the entire attachment line of the models over 
the freestream Reynolds number range tested (consistent with Poll’s fourth regime). 
When all end disturbance sources are removed, transition occurs on the attachment 
line independently of freestream noise levels and in agreement with the first regime. 
Finally, with the addition of small roughness elements on the cylinder attachment 
line, transition occurs at lower values of the Reynolds number, depending on both 
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the roughness height and the Wind-tunnel noise level (characteristic of the second 
and third regime). 


The evaluation of /(, 




/dx) x ^ 


we have 


^Adujdx)^ 


II = R = Vt - ( 27 ) 

V "♦ \/i'*(du e /dx) x= 0 

Introducing the Reynolds number based on the freestream conditions and the local 
diameter D of the leading edge, we have : 

*• = 7 1 CM) 

Defining now the freestream velocity in terms of its normal and tangential components 
to the leading edge we write : 

* u oc,n ,/ u oo,n /")n\ 

cos A = — — or Voo = V&) 

Voo cos A 

Since = v 2 = V e = sin A we also have 


— sin A 


Replacing equations (29) and (30) in (28) we obtain : 


— sin A - . \/ Roo,d 

v * -/cosA ^{dujdx) x ^ 
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We still need to define the pressure gradient in terms of known quantities. We 
know that in the stagnation region we can write 8 



Since 

^ oo,n 

^ OO , fl 1 ») rfi 

7 Hgas * oo 

we can replace it in the previous equation to obtain : 



According to the modified Newtonian theory, we can write 


( 33 ) 

( 34 ) 

( 35 ) 

( 36 ) 

( 37 ) 


Pc - Poo = (PU2 - Poo) cos2 A 


( 38 ) 


where 


^ = /(M oo,7) 
Poo 


( 39 ) 


is given by the Rayleigh Pitot tube formula. We still need to evaluate the temperature 


ratio : 



Poo _ Pi _ /;(Af ooCosA,7) (normal shock relation) 

P'2 P'2 


( 41 ) 
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To T 2 

~~ = — = /;/(Moo cos A, 7 ) (normal shock relation) (42) 

too i\ 

A last step needs to be completed in order to evaluate the kinematic viscosity ratio : 

(43) 


I r'oo I Poo P* 

V V * \ P* PoO 


Using Sutherland’s law of gas viscosity, we can calculate the ratio Poolp* as 

- = (W££ aM 5=110 K 

P* \ * * / 


(44) 


In order to calculate p,/ poo we will assume a constant pressure through the boundary 
layer. Therefore, the density ratio can be expressed in terms of the temperature and 
pressure ratios. 


P* P* ^00 Pe Too 

Poo Poo Poo T t 


(45) 


Since we know how to evaluate the pressure ratio, the only quantity which remains 
to be evaluated is the reference temperature. According to Poll 13 , we have : 


T t = T e + 0.10(7;, - 7 ;) + 0.70(7; - T e ) 


(46) 


with the definition of the recovery factor 

Tr-T' 


r = 


and Tq = Too [ 1 + 


(' + 2 t 1m ~) 


7’o - T e — ~ ’ 2 

We have now identified all the variables to evaluate R ,. 
Finally, we can rewrite R, in the following way : 

fMZ* f** sin A 


(47) 


(48) 


(ifcO-*))' 

It appears immediately that increasing the diameter or the unit Reynolds number will 
be unfavorable for the transition at the attachment line, since R, will increase. 


20 


An increase in sweep angle is also unfavorable since both p e /poo <md Tg/Toc 
will decrease and thus increase the value of /?.. 

The effect of the wall temperature can also be analyzed through the viscosity ratio. 
When the wall temperature is increased, the reference temperature will increase and 
both the density ratio and the dynamic viscosity will decrease. Increasing the wall 
temperature along the attachment line will have a favorable effect. 


The evaluation of r/ 

As defined by Poll, we have : 


= ( * — ) 

\ {du € /d^) Xst0 ) 


1/2 


Making use of the previous derivation of ( du e /dx) x=Q we have 

/ \ ,/2 
| n^y/W^T, 
v = 


V 


Pr 




where 


11 . = 

Pref \Tref) T f 


T n f+S 


p Tt f \ Jr. ( l If + S 
with S=110 K, T ref = 273.1 K and /i ref =0.1716 l(T^kg/(m s). 


(49) 


(50) 


(51) 


Transition due to cross flow instability 

Let’s assume that transition along the attachment line has not taken place and 
the flow is still laminar. In the region close to the attachment line, the streamlines 
at the edge of the boundary layer are highly curved (fig-1 1). Along any of these 
streamlines there is a favorable pressure gradient and we should expect that the flow 
will be stable to small disturbances. However, it turns out that this is not the case. 
Poll 12 explains that by simply sweeping the cylinder shown in figure 1 1, it is possible 
to destabilize the boundary layer since the sweep leads directly to the formation of 
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velocity profiles presenting inflection points (fig.15 and fig.16). To properly determine 
the onset of transition due to cross-flow instability, the most reliable criteria seems 
to be the characteristic local Reynolds number X. 

X is based on the cross-flow velocity profile shown in figure 16 and corresponding 
to the case e = 8 - f. X is defined as : 

V 

X _ ^mar /v Q OiC mnl (52) 

// 

Where C maI is the maximum cross-flow velocity component in the profile corre- 
sponding to t = 0 - | and Z 00i c mat is the point where this velocity has dropped to 
1% of its maximum value. An excellent correlation seems to exist between the criti- 
cal value of X and the streamwise shape factor //„ as shown in figure 17 (excluding 
situations where the streamwise flow is unstable to small disturbances). Unfortu- 
nately, this kind of transition criteria requires a detailed boundary layer calculation 
in the leading edge region that the waverider code does not provide. Therefore, it 
was decided to take into account the destabilizing effect due to sweep by simpler 
means. In 1985, Bowcutt 2 curve fitted experimental data obtained on blunt flat plate 
with a swept leading edge in such a way that the transition Reynolds number in the 
streamwise direction for a swept leading edge (Rc xt )\ could be related to the unswept 
value (/?e Z f) A=u by a function of the sweep angle only (fig.18). Since {Re xt ) A=0 will 
be determined as explained later on and the sweep angle distribution is known, we 
are able to calculate ( Re x t) A for each streamline. Of course this procedure will be 
applied only for the streamlines where it was found that the flow was still laminar 
at the attachment line. When the flow was found to be turbulent at the attachment 
line, this criteria does not apply since the flow will remain turbulent on the streamline 
from the leading edge to the trailing edge of the vehicle. 

The possibility of relaminarization 

The reversion of a turbulent boundary layer to the laminar state in a region 
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of strong favorable pressure gradient is a well established phenomenon in two- 
dimensional and axi-symmetric flows, _ 

In a simplified analysis, Poll 11 demonstrates that the conventional relaminarization 
parameter reaches its peak value very close to the attachment line. Therefore, 
turbulence resulting from either attachment line contamination or cross-flow instability 
is unlikely to be relaminarized in the strong favorable pressure gradients which exists 
in the streamwise direction near the leading edge of a swept wing. 

Although relaminarization is not likely to occur in the streamwise direction, 
we will underline here that in the spanwise direction (along the attachment line) 
relaminarization might occur since the attachment line cannot support bursts of 
turbulence when R , drops below 245 (cf. section 2 topic 1). 

Transition in the streamwise direction 

In the previous work of Bowcutt, Corda and McLaughlin a correlation based 
on transition Reynolds number of a large bulk of measurements for sharp cones 
(fig. 19) was used. This correlation, when compared to a larger set of experiments 
(fig. 20) grouping flight measurements and quiet wind tunnel experiments, shows a 
good agreement with the wind tunnel data correlations presented in figure 20. 

Therefore, due to the well known adverse effect of wind tunnel radiated noise on 
the transition Reynolds number, it was decided to try another correlation based on 
ballistic tests performed on sharp cones by Sheetz 15 . These data can be correlated 
in such a way that the transition Reynolds number could be expressed as a function 
of the edge Mach number : 

Rcjt = 640o(A/e) 36f ' (53) 

This relation, also presented on figure 19, shows in the range 6 < Me < 16 a better 
agreement with flight data and the c ^ method than Bowcutt s curve fit. 


23 


Both correlations do a very poor job of predicting transition when the edge Mach 
number drops below 5. As stated by Malik 16 , at low hypersonic speed ( M < 7) 
the first Tollmien-Schlichting mode dominates for adiabatic wall conditions. For cold 
walls, the role of the second mode becomes more important due to the destabilizing 
effect of cooling. We will see later on that when radiative equilibrium is obtained at 
low hypersonic Mach number (A/oo= 5), we can always consider that T W /T Q m o.s 
is a good estimate. Therefore, for Moo < 7 and T w /T 0 * o 8» we will consider a 
constant transition Reynolds number of 15 million. 

At higher Mach number 10 < M^ < 20, T w /T 0 of the order of 0.1 or less 
is a good estimate for radiative equilibrium wall temperature. For these conditions, 
we know that the second mode is dominant and Sheetz’s transition criteria does a 
relatively good job for M e > 6. 

In fact, for 6 < M e < 16 , Sheetz’s and Bowcutt’s criteria can be seen as the 
upper and lower band bracketing almost all the flight data. 

Section 3 The heat transfer evaluation 

The heat transfer evaluation on the complete vehicle will be separated in three 
distinctive regions : the nose region, the attachment line along the leading edge, and 
the upper and lower surfaces. The heat transfer rates in the two first regions will 
be obtained through analytical correlations developed by Tauber 17 . These relations 
apply in the flight regime where boundary layer theory is valid and the flowfield is 
in equilibrium. The peak hypersonic aerodynamic heating experienced by missiles or 
manned vehicles always occurs in the continuum flow regime. However, for small 
nose radius or wing leading edge radii, these regions of the vehicle can experience low 
density flow phenomena which alter the heating. One parameter that characterizes 
low-density flows is the Knudsen number, Kn, which is the ratio of molecular mean 
free path to a characteristic length. The increase in stagnation point heat transfer 
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coefficient with Knudsen number was calculated by Moss and is shown in figure 21. 
We see in this figure that there is a good agreement between the direct Monte-Carlo 
simulation, the viscous shock layer solution with no velocity slip and the viscous 
shock layer solution with velocity slip for Kn less than 0.05. Beyond that point, 
the viscous shock layer solution ceases to be valid and overpredicts the heat transfer 
calculated by the DMSC. Since the correlations developed by Tauber will assume 
continuum flow, we need to ensure that the Kn number will always be less than 
0.05. If a characteristic length R min , equal to the minimum leading edge radius is 


assumed, then ; 


/\ n — = 0.05 OT Rmtn — 20AQO 

fimin 


Considering the definition of the mean free path, 


ArjO — 


with Hoc ~ 


Poo 


\/2ird 2 n 0 o kToo 

where d = 3.789 10~ 10 m for pure N 2 and k = 1.380 54 1 0 - 23 J/K. 


(54) 


(55) 


R min can be rewritten as : 


R 


min 


20 fcT'oo 

V^ird 2 Poo 


(56) 


Rmin represents the smallest radius that we can consider at a given altitude in order 
to avoid low density effects. Therefore, we will always make sure that the smallest 
radius on the leading edge is greater than this minimum value. When this condition 
is not met, the entire waverider will be rejected by the optimization procedure. 

The heat transfer rates on the upper and lower surfaces will be evaluated assuming 
a 2D flowfield along the streamlines through the reference temperature method and 
the Reynolds analogy. 

Finally, we will consider that the surface of the vehicle is a fully catalytic wall 
where no boundary layer mass addition is present. Assuming a fully catalytic wall is 
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a conservative approach and makes the heat transfer rates independent of the choice 
of heat shielding material. 

Convective heat transfer evaluation in the nose region 

Tauber 17 demonstrated that the heat transfer rate at the stagnation point of an 
axisymmetric body can be expressed as : 

<«»>•" = ('“'»-) V “(l-^)(7|) (Watt/m 2 ) (57) 


• Ry in meters 

• p in kg/m 3 

• V in m/s 

This relation is derived from a similar solution of the boundary layer equations 
where air is assumed to be a binary mixture of “air atoms’ and “air molecules” 
in thermodynamic equilibrium. The equations are then numerically integrated using 
the composition and thermodynamic properties of air in chemical equilibrium and a 
unit Lewis number. 

The resulting correlation, almost identical to the Marvin and Deiwert correlation, 
is compared with experimental data in figure 22. Up to a speed of 14 km/s all four 
correlations lie within the data spread, which is about ± 25%. 

By using the relation given by Hamilton 18 , we can obtain the 3D stagnation heat 
transfer from the axisymmetric value. 

(Mw = «.).„• ( 58 ) 

K = for R y < R x (59) 
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( 60 ) 


K - ~ for R z < R y 

tty 

Note that K ranges always between 0 and 1. For K=1 we have an axisymmetric nose 
and for K=0 we have a cylinder. R x has to be seen in our case as the radius joining 
the stagnation point and the second point of the planform (fig.7). 

The wall enthalpy was calculated assuming a calorically perfect gas at the wall. 

h w = Cp T w ( 6 1 ) 

Since we are considering a fully catalytic wall, a high pressure region and wall 
temperature not exceeding 1900 K for any realistic configuration, we see that the 
assumption of a calorically perfect gas at the wall is reasonable. Only the effect of 
vibrational excitation is neglected in that case. 

When the stagnation heat transfer at the nose is known, the heat transfer at the 
second point of the planform is calculated according to the value of K. 

• For 0 < K < 0.5, the heat transfer distribution between the stagnation point and 
the second point will be close to the heat transfer distribution existing on an 
unswept cylinder of radius R x . For this range of K, we can write 

= 0.7 (cos A, =2 ) 3/ ' 2 + 0.3 (62) 

{*iu>) 3 tag 

where Aj =2 represents the sweep angle at the second point. 

• If 0.5 < K < 1.0, the heat transfer distribution will be closer to the heat transfer 
distribution of an axisymmetric body and we can use 8 : 

iisW = cos A,., <«> 

{Qw)atag 


27 


Convective heat transfer evaluation at the leading edge 

The heat transfer along a cylindrical leading edge of a finite length wing can 
be approximately calculated using the expression developed by Rubesin and given 
in reference 17. 


(</«-)/«. = ((4w)c y i + sin 2 A (4w)j p ) / cos 


a 


(64) 


Since we will only deal with a waverider at zero angle of attack and the attachment 
line of the leading edge, we can set a to zero. But we must realize that this equation, 
for highly swept leading edges, is limited to small angles of attack since the stagnation 
line moves from the cylindrical leading edge onto the lower surface of the wing as 
the angle of attack increases. 

The heat transfer rate for the infinite swept cylinder is calculated with 17 : 

1/2 


= *(•-£)(*) 

* cosA(l -o.lSsin 2 A) ( Watt/m 2 ) 


(65) 


V 2 

haw — hoo H — — (i-o.i8sin 2 A) 


For the inclined laminar flat plate we have 17 : 


(4w)f p = ( 2.42 to -5 ) V£ 2 1 


hw \ f Poo cos A \ 

h-aw ) \ ‘S / 


1/2 


sin 6 ( Watt/m 2 ) 


( 66 ) 


(67) 


with 6 = — — A (68) 

2 

V 2 

and h aw = ^oo + -^(l-o.issin 2 A) (69) 

2 ' 

As in the case of the axisymmetric stagnation point correlations, these two relations 
are derived from exact similar solutions and have assumed the same air properties. 
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To make sum that the heat transfer predicted at the second point of the plantain 
by the swept cylinder correlation is equal to the one predicted by the stagnation 
point correlation, an artificial origin is created for the laminar flat plate correlation. 
Therefore, no discontinuity in the heat transfer at the second point of the planform 
shape will exist. Note also that this manipulation is somewhat conservative from the 
heat transfer point of view since it will slightly raise the laminar flat plate contribution 
of the downstream points. 

When transition along the leading edge is detected, the evaluation of the flat 
plate contribution has to be done through a turbulent correlation. This correlation, 
developed by Tauber 17 is obtained by using : 

• the Reynolds analogy, 


. the Blasius incompressible, turbulent, skin-friction equation for Reynolds number 
less than 10 million, 

Cj_ _ 0-0296 (7i) 

2 ~ (Re,? 2 


. the Schultz-Grunov equation for Reynolds number above 10 million, 

Cj = 0-185 ( 72 ) 

2 (hgwRe.f™ 


The denotes that the quantity is evaluated at the conditions corresponding to the 
reference enthalpy h.. 


ht — 0.18/l a „i + 0.32 h e + 0.5h w 


(73) 
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Using the thermodynamic and transport properties of air in chemical equilibrium, 
Tauber calculated the heat transfer rate at altitudes from 3 to 60 km and speed from 1 
to 1 1 km/s. His correlation which is accurate to ± 15%, may be expressed as follows : 


• for 1500 < V iyo < 3960 m/s 

{iw) fp = ( 3 . 7210 -*) Vi” 



hw\ (poo sin 2 6 cos 2 22 <5)° 8 
W (3-W«) 0 ' 2 (^/555) 0 ' 25 

(Wattjm 2 } 


(74) 


and for Voo - 3960 m/s. 

. 37 / h w \ (poo sin 2 S cos 2 62 <5) 0 8 

(Watt/m 2 ) 


(75) 


with fi = - - A (76) 

2 

V 2 

and h aw = h ^ + -^(l-o.ll sin 2 A) (77) 

When transitional heat transfer values have to be evaluated, we will assume that the 
flat plate contribution will vary linearly with the distance s between the laminar value 
at the onset of transition and the turbulent value at the end of transition. 

When the flow undergoes reverse transition or relaminarization, the same hypoth- 
esis is made but the heat transfer rate will vary from a turbulent value at the beginning 
of relaminarization to a fully laminar value at the end of relaminarization. In order to 
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calculate {q w )/ p after relaminarization, a new origin has to be defined for the laminar 
correlation. This is done by imposing that the flat plate heat transfer contribution 
at the beginning of relaminarization be the same for the laminar and the turbulent 
correlation. Since the relaminarization length is predicted by our transition criteria, 
we are able to evaluate the laminar heat transfer value at the end of relaminarization 
and all the other downstream points. 

Convective heat transfer evaluation on the waverider surface 

In order to evaluate the viscous forces on the waverider, both Corda 3 and 
McLaughlin 4 have used the so called reference temperature method to predict the 
skin friction coefficient of laminar and turbulent flows. The equations used in their 
work are : 

• For laminar flow 

1 / 7’* \ l 1 " -1 )/ 2 

<7,.„ = 0.664(Ke,)-^— J 0 8) 


with 


Re, = 


PeVeS 

Pe 


(79) 


T* is defined by 


( — \ = 1 + 0.032 + 0.58 ( ^ - l) (80) 

\7oc / \ Joo / 


and u?=0.75 which is the exponent of an assumed exponential variation of the 


viscosity, where // is defined as 



(81) 


• For turbulent flow 


0.0592 

C,, "‘ ~ (Re,.? 2 
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(82) 


1 



where 


Re st = tYs2i 

P. 


( 83 ) 


The present study will also use the reference temperature method. The Stanton 
numbers are calculated through the Reynolds analogy. The Stanton number is defined 
as follows : 


St 




PeVe(haw ^tt>) 

If a calorically perfect gas is assumed, equation 84 can be written as : 

(iu,) 


(84) 


St = 


PeVe c p (T aw T w ) 


(85) 


Making use of the definition of the recovery factor, we obtain the laminar or the 
turbulent heat transfer rate through : 

(<?u>) = t p e V c St Cp (Tq Tiw) (86) 

where r = y/~Pr for laminar flow and r = Pr ^ for turbulent flow. 

Section 4 Evaluation of the equilibrium radiative heat transfer 

We have seen in the previous chapter how to calculate the heat transfer brought 
by convection to the surface of a hypersonic vehicle. In steady or quasi-steady flight, 
pan or the entire convective heat transfer brought by the fluid will be radiated by 
the vehicle surface. Therefore, heat transfer due to radiation constitutes an important 
natural cooling phenomenon that must be taken into account. 

The radiative heat transfer emitted by a surface which has a wall temperature 
T w is expressed as : 

(9w)rad = t ° PJ — 72,) (87) 
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with 


. e = surface emissivity 
. a - Stefan-Boltzmann constant 

For a given wall temperature of the surface the heat transfer balance at the wall is : 

(9u'(r u p )) balance = (Vw(T w )) c <n ivtci ~ (4u>(^))r<id 

To maintain this given wall temperature, the heat transfer balance has to be 
absorbed either by an active cooling device or by the vehicle structure. Therefore, 
the hea, transfer balance represents the amount of cooling necessary a. a given point 
of the vehicle to maintain an imposed wall temperature. 

If the vehicle structure doesn’t provide any heat sink ((9u-)b a i an ce _ °)’ 
temperature will reach a value which is such that the heat transfer losses by radiation 
balance the heat transfer due to convection. This temperature is called the radiative 
equilibrium wall temperature and is calculated using the following relation . 

= = 0 (89) 

In order to find this radiative equilibrium wall temperature, a aero finding 
algorithm due to Brent 19 was used since the convective and radiative heat transfer 
are both non-linear functions of the wall temperature. In other terms, it is necessary 
,o iterate at each point of the waverider surface to obtain the radiative equilibrium 
wall temperature and the corresponding convective heat transfer. 

The waverider code was set up in such a way that the maximum wall temperature 
at the leading edge and on the upper and lower surface could be independently 
specified by the user. The surface emissivity was in all cases arbitrary set to a value 

of 0.8. 
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When the radiative equilibrium temperature exceeds the maximum wall temper- 
ature allowed by the thermal protection system, the wall temperature is set to this 
maximum value and the convective heat transfer is recalculated using that value. 

Therefore, the user can rapidly evaluate the amount of cooling necessary to 
maintain its maximum wall temperature by substracting the radiative heat transfer. 
The location and the extend of the surfaces which have to be cooled can easily be 
visualized by temperature contours on the waverider surface. 


Section 5 Viscous forces integration 

The procedure to integrate the viscous shear stresses is identical to what was 
used and described in detail in reference 4. However, differences in the transition 
front will appear since Bowcutt 2 , Corda 3 and McLaughlin 4 have considered a simpler 
transition model than the one used in this study. In particular, they have assumed a 
constant and identical wall temperature on the upper and lower surface to calculate the 
shear stress distribution. In this study, the shear stress at each point of the waverider 
surface is evaluated at the corresponding radiative equilibrium wall temperature or 
the maximum wall temperature. The code will therefore provide a more accurate 
evaluation of the viscous forces. This will be particularly true at high Mach number 
where the difference between the temperature of the upper and lower surfaces will 
tend to increase. 


Section 6 Definition of a thermal protection system 
Leading edge 

Three different types of generic leading edge concepts are listed below for an 
application such as the NASP. 
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• radiative or "passive” concept : the cooling relies simply on radiative cooling. 

It is a reliable and fail-safe solution which doesn’t require coolant or pumps and 

controls. 

. heat-pipes or “semi-passive” concept : the cooling relies on radiative cooling 
and heat pipes which transport heat to a cooler region of the wing where it is 
dissipated radiatively. 

• heat-pipes and active cooling or “semi-active” concept : the semi-active concept 
uses heat-pipes with active cooling. The active cooling is used during the ascent 
and cruise portion of the flight. During the descent, since the heat transfer will 
be much lower, active cooling is not necessary and the concept is used in a 
semi-passive way. Since this concept requires some form of controls, it has to 
be considered as less reliable than the two first solutions. Nevertheless, it might 
be the only solution able to support high heating rates with no damages. 

It should be also noted that pure active cooling, such as channel flow impingement 
or transpiration, could be used but offers even less reliability and requires heavier 
control systems. Therefore, these systems have to be disregarded for the leading edge. 

It has thus been decided that the most suitable system for our waverider should be 
the semi-active concept. In reference 20, a carbon-carbon / refractory metal heat-pipe 
wing leading edge is studied in the context of the NASP. 

This concept is designed to be redundant and fail-safe and to operate at temper- 
ature as high as 1900 K. This concept uses internal radiative cooling to supplement 
external radiative cooling during the ascent portion of the trajectory. 

The design, shown in figure 23, features a carbon-carbon primary structure with 
refractory metal chordwise and spanwise. The heat-pipes are made of tungsten with 
lithium used as the working fluid. The heat-pipes are completely embedded within 
the carbon-carbon structure. The spanwise heat-pipes are much smaller in size than 
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the chordwise heat-pipes and are located only in the stagnation region. 

Although carbon-carbon retains strength up to 3000 K, the maximum reuse 
temperature is limited to approximately to 1900 K by the coating required for anti- 
oxydation at elevated temperature. Note that, due to its properties, the carbon-carbon 
also offers some ablative protection in the event of a heat-pipe failure. The maximum 
temperature as a function of the stagnation heat transfer is presented in figure 24 for 
3 cases : no heat-pipes, no internal cooling, and heat-pipes with internal cooling. 
For the first case, temperatures were approximated by the radiation equilibrium 
temperatures. Only the third case is able to maintain the maximum temperature 
below the temperature limit over the predicted heat transfer range for the ascent 
trajectory. As shown in the figures, the thermal limitations of the system are a 
maximum temperature limit of 1900 K and a maximum heat transferrate of the order 
of 10.3 MW/m 2 (900 BtuJ ft 2 5 ). Due to its relative complexity, a minimum value 
of the leading edge radius must be imposed and was set according to the study 20 to 
1.27 cm (0.5 inch). 

Waverider surface 

As in the case of the leading edge, different choices can be made according to 
the maximum temperature, the density, the reliability, etc., of the thermal protection 
system. Basically, the choice reduces to the three following categories of material 1 : 

• carbon-carbon composites (1650 < T max < 2200 C) 

• advantages : 

• low specific gravity. 

• high modulus. 

• low thermal conductivity. 

• strength increases with temperature. 
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• excellent thermal and mechanical shock resistance and toughness. 

• disadvantages : 

• oxidation is a problem but rates arc not catastrophic in oxidizing 
atmosphere. 

• due to corrosion, coating is required for long life performance and 
the maximum temperature falls to 1650 C due to the coating. 

• ceramic composites (1400 < T fnai < 1600 C ) 

• advantages : 

• excellent oxidation resistance. 

• specific gravity, specific strength and specific modulus values are excel- 
lent. 

• disadvantages : 

• thermal stresses and thermal and mechanical shocks are more severe than 
with the carbon-carbon composites. 

• superalloys of titanium ( T max < 1100 C ) 
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Chapter 4 Results and discussion 

Section 1 Choice of flight trajectories 

In order to determine test cases of interest corresponding to an application similar 
to the NASP, it was decided to choose two constant dynamic pressure trajectories. 
Values of 0.2 and 1.0 atm were chosen. The low value corresponds to what is believed 
by the airframe community to be the maximum admissible from a structural point of 
view. The high value corresponds to what is desired by the engine community to 
obtain the best performance of the scramjet. For both dynamic pressure trajectories, 
the waverider performances are analyzed at four different Mach numbers (5, 10, 15, 
20). The flight conditions corresponding to each dynamic pressure and each Mach 
number are given in Table 1. Pressure and density were extracted from reference 21, 
in order to match the desired dynamic pressure as close as possible. The velocity 
ranges from 1500 m/s to 6800 m/s and the unit Reynolds number varies from 0.35 to 
9.35 million per meter. The low dynamic pressure flight path is situated between 30 
and 50 km of altitude and the high dynamic pressure trajectory is situated between 
20 and 40 km. The two flight path are showed on a Mach-Altitude map on figure 25 
which also features the “official “ NASP trajectory. 

Section 2 Choice of other parameters 

The choice of an altitude and a Mach number does not represent the entire set 
of parameters which have to be chosen to run a test case. Another set of important 
parameters have to be defined in the input data file and we will review them one 
by one. 

• The objective function for the optimization : a choice can be made between 
L/D, Cd, Volume or Volume efficiency. In this study, we will only consider 
waveriders maximized for L/D or Volume. 
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• The cone angle : associated with a given Mach number for the calorically perfect 
gas or with given freestream conditions in the case of a gas in chemical and 
thermodynamical equilibrium, fixes the streamlines of the generating conical 
flowfield. The choice of the cone angle has an impact on the optimized objective 
function. It was demonstrated 2,3 that for each flight condition there exists a 
certain value of the cone angle leading to the best optimum for the optimized 
objective function. Nevertheless, in the present study, the notion of best optimum 
was disregarded and an arbitrary well chosen value of 7.5 degree was considered 

in all cases. 

• Type of gas : the program can be run with a generating flowfield corresponding 
to a calorically perfect gas or a gas in chemical and thermodynamic equilibrium. 
Since it has already been demonstrated 4 that the difference between the two 
cases is small for low altitude and Mach number smaller than 25, this study has 
considered a calorically perfect gas for the generating flowfield in all cases. 

• Length of the vehicle : all test cases have considered a 60 meter long waverider. 

• Geometric constraints : were placed on the vehicle maximum box size (0.4) and 
minimum slenderness ratio (0.075). 

• Transition at the attachment line of the vehicle : the transition criteria at the 
attachment line depends on the surface roughness of the leading edge. For all 
the test cases we have considered a “rough” leading edge. This choice will be 
justified later on. 

• Maximum wall temperature : according to chapter 3 section 6, a maximum wall 
temperature of 1900 K was chosen for the leading edge and the waverider body 
although they might be covered by different materials. 

• Cooling and minimum radius : for the two test cases at Mach = 5, cooling is 
not necessary if carbon-carbon is used to protect the leading edge. Therefore, 
cooling was not included for in these two cases and the minimum radius of the 
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leading edge was set to 1 mm. For all the other test cases, cooling is required 
to maintain a maximum wall temperature of 1900 K. According to the geometry 
of the active cooling device, the minimum radius was set to 10 mm. Although 
the choice of cooling can be made in the input data file, the constraint on the 
minimum value of the radius has to be changed in the code itself. 

The following options have to be modified in the code : 

• Waverider’s carving : the amount of carving is set in the program according to 
equation (1). For all cases, the parameters used to define the carving function 
were set to the values given in chapter 2 section 3. 

• Transition in the streamwise direction : for < 7 the critical Reynolds number 
was set to 15 million. For > 7 the critical Reynolds number was calculated 
according to Sheetz’s criteria. The user may choose to go back to Bowcutt’s 
criteria if necessary. 

• Surface emissivity : the surface emissivity of the waverider was assumed to be 
constant with temperature and set to a value of 0.8 for all surfaces. 

The complete definition of each test case is summarized in Table 2 and 3. 

Section 3 Discussion of the results 

Waverider geometry 

The waverider geometries are presented in figure 26 for the low dynamic pressure 
and in figure 27 for the high dynamic pressure. Due to the high values of the Reynolds 
number, both Mach 5 waveriders are similar to the caret wing which represents the 
optimum inviscid waverider. For a constant dynamic pressure trajectory, as the Mach 
number increases the Reynolds number decreases. Therefore, at higher Mach number, 
the geometry is less likely to resemble a caret wing. Note also, that the geometries 
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for Mach greater than 5 do not change very much; there is a reasonable amount of 
Mach number independence at high Mach number. 

The principal dimensions are given in Table 4. As expected, for increasing Mach 
number and constant cone angle, the shock angle as well as the span will decrease. 
It is interesting to note that on both trajectories, the base height increases with Mach 
number. The effect of Mach number on the span and the base height results in an 
increase in volume efficiency with Mach number along a constant dynamic pressure 
trajectory. For Mach number less or equal than 10, the effect of altitude on the 
geometry is nonexistent : the two shapes are almost identical. For Mach greater or 
equal to 15, the effect of altitude on the geometry is small and acting only on the 
span of the vehicle. 

Leading edge results 

For a given geometry and a given carving law we can define a radius and a 
sweep angle distribution along the leading edge. These quantities, when combined 
with the freestream properties give us the ability to calculate the transition criteria R t 
and the boundary layer characteristic length rj. These four quantities are presented 
on figures 28a, 28b, 28c, 28d for the Mach 5, low dynamic pressure case. Figure 
28e represents a flag function of the flow nature along the leading edge. When the 
flag is set to zero, the flow is laminar or has relaminarized. When the flag is set to 

1, the flow is in a transitional or relaminarizing state. Finally, when the flag equals 

2, the flow is fully turbulent. According to the values of R , and rj in this case, 
we see that the flow is purely laminar all along the attachment line of the leading 
edge. The next picture (fig.28f) describes the variation of the pressure gradient along 
the leading edge. This pressure gradient is calculated according to the modified 
Newtonian theory. Finally, figures 28g and 28h show the convective heat transfer 
(in megawatts/m 2 ) and the radiative equilibrium wall temperature distribution (in K) 
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along the leading edge. Figure 29 illustrates the same features for the high dynamic 
pressure case. As expected, due to a factor 5.6 increase in unit Reynolds number, 
R, increases by a factor 2.4. The increase in unit Reynolds number produces the 
opposite effect on r/ which drops by more than a factor 2. As a result, we see that 
for almost identical radius and sweep distribution, the low dynamic pressure case 
presents a fully laminar leading edge and the high dynamic pressure case a fully 
turbulent one. This major difference leads to an increase of the radiative equilibrium 
wall temperature which is of the order of 50 K. Since we have chosen a leading 
edge made of carbon-carbon with a maximum temperature of 1900 K, the choice 
for a non active cooling option was justified for both cases. Looking to the values 
of 7) we can also justify the choice of a “rough” leading edge surface. For the low 
dynamic pressure, t] is of the order of 0.15 mm and for the high dynamic pressure 
7) is even smaller (0.06mm). We have previously defined our rough leading edge 
by the ratio d A/ > 2.4. Therefore, at low/high dynamic pressure any obstacle larger 
than 0.36/0.14 mm along the attachment line would represent for the flow a major 
disturbance. These values justify a posteriori the choice of a “rough” leading edge. 

If we now increase the Mach number to 10 (fig.30 and 31), the combination 
of larger radius along the leading edge and smaller unit Reynolds number produce 
R , values which are slightly lower than for the Mach 5 case. In this case, we find 
the same situation at Mach 10 and at Mach 5. The flow is laminar/turbulent at 
low/high dynamic pressure. Due to the increase in the velocity, the heat transfer has 
substantially increased in both cases. At low dynamic pressure, the wall temperature 
now ranges from 1400 K to 1900 K. Due to the turbulent nature of the flow, at high 
dynamic pressure, the wall temperature rises a little bit and varies between 1600 K 
and 1900 K. Note that in both cases, active cooling becomes necessary but can be 
limited in a small region near the nose. The critical roughness size is now of the 
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order of 0.96/0.43 mm for the low/high dynamic case. The assumption of a rough 
leading edge is still reasonable at Mach 10. 

At Mach 15 (fig.32 and 33), only the low dynamic case presents a fully laminar 
attachment line. At this point, the high dynamic pressure case requires active cooling 
all along the leading edge. The low dynamic pressure requires active cooling only 
in a small region near the nose. The average value of r) continues to rise for both 
cases when compared to the Mach 5 and 10 cases. Here, with a critical roughness 
height of 1. 6/0.8 mm for low/high dynamic pressure our assumption might not be 
totally valid anymore. Note that, at Mach 15 and 20, this is of no consequence at 
all. Assuming a smooth leading edge will only raise the critical value of needed 
to provoke transition. The rough leading edge option has therefore to be seen as a 

worst case scenario. 

Finally, at Mach 20, both leading edges are fully laminar and require active 
cooling all along the leading edge. At high dynamic pressure, the convective heat 
transfer is maximum at the nose and of the order of 9.5 MW/m 2 . The radiative 
heat transfer at 1900 K is of the order of 0.6 MW/m 2 . Therefore, the active cooling 
system should be able to absorb 8.9 MW/m 2 . This value, for the chosen active cooling 
device, is well within the maximum heat transfer of 10.3 MW/m 2 admitted by the 
system. For the rest of the leading edge, the convective heat transfer is much lower 

and ranges between 1.0 and 2.0 MW/m . 

Transition front and heat transfer characteristics on 
the waverider surfaces. 

Once the flow nature is known along the attachment line, we have in fact our 
“initial flow properties” at the first point of each upper and lower streamlines. When 
the flow is turbulent along the attachment line we know that it will stay m that 
condition all along the streamline up to the trailing edge of the vehicle. When the 
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flow is laminar at the leading edge, it might undergo transition along the streamlines 
according to the value of the local Reynolds number. 

The transition front on the upper and lower surface of the vehicle are given for 
our eight cases in figures 36 to 43. On each figure, the black surfaces are indicative of 
laminar region. When the flow is turbulent, the details of the grid are apparent. Due 
to their different way of construction, the upper and lower surfaces are different. The 
upper surface is always represented by a matrix of 51x51 points (some overlapping 
on the planform curve). On the other hand, the lower surface is represented by 51xK 
points, where X may be any number between 2 and 100. The value of K, which 
represents in fact the number of stations in the flow direction, varies according to the 
choice made by the optimization algorithm on the five initial points (see chapter 2 
section 2 or reference 4 for more details). The lines oriented in the flow direction 
on the upper surface or with a small angle on the lower surface are the streamlines 
of the flow on both surfaces. 

At Mach 5, and low dynamic pressure (fig.36) we know that the flow is laminar 
all along the leading edge. Meanwhile, the unit Reynolds number is so high that 
almost the entire surfaces are wetted by turbulent flow. At high dynamic pressure 
(fig. 37), the flow is fully turbulent at the leading edge and as a consequence both 
surfaces are wetted by turbulent flow. 

At Mach 10 and low dynamic pressure (fig.38) we can see that the laminar 
flow region is more important and representative of the combined favorable effect of 
decreasing the unit Reynolds number and increasing the Mach number. At high 
dynamic pressure (fig.39), the flow is still turbulent along the leading edge. In 
consequence, both surfaces are still wetted by turbulent flow except near the nose. 

At Mach 15 and 20 (fig.40 to 43), the leading edges are laminar in all cases. 
Therefore, at a given Mach number we can compare the effect on the transition front 
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due to a change in unit Reynolds number. Or, at a given dynamic pressure, we can 
visualize again the combined favorable effect of decreasing the unit Reynolds number 
and increasing the Mach number. Note that at high dynamic pressure (fig.41 and 43), 
most of the lower surface is still wetted by turbulent flow. It is interesting to note 
that for the high dynamic pressure trajectory, most of the flow on the lower surface 
is turbulent from Mach 5 to Mach 20 (fig. 37, 39, 41, 43). Now that the flow nature 
has been determined along each of the upper and lower streamlines we can evaluate 
the laminar and turbulent heat transfer characteristics along these lines. Since the 
high dynamic pressure trajectory always presents higher heat transfer we will focus 
on that particular one. 

For each Mach number, the convective heat transfer on the upper and lower 
surface is presented along different streamlines (fig.44 to 47). The main purpose of 
these graphs is to give to the reader an order of magnitude of the convective heat 
transfer involved at each Mach number. Note also that the curves are numbered from 
1 to 51. The curve number 1 is the centerline streamline and the curve number 51, 
represented by one point, is the wing tip. By joining the first point of each curve, 
we obtain the heat transfer distribution along the attachment line of the leading edge 
presented earlier. By choosing a given streamline, we can follow the evolution of the 
heat transfer as the flow travels from the leading edge to the trailing edge. The Mach 
15 case is rather interesting (fig.46) since it illustrates fairly well on the upper surface 
the effect of sweep angle on the critical Reynolds number as explained in chapter 3, 
section 2, topic 2. Since the heat transfer is always higher on the lower surface, the 
next set of figures (48 to 51) presents the wall temperature distribution only on the 
lower surface. To allow more details, the span has been artificially expanded here. 
The main purpose of these graphs is to determine what kind of thermal protection is 
needed and if eventually active cooling would be required in certain zones. 
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Obviously, at Mach 5 (fig.48), most of the surface is at a wall temperature of 
800-900 K. 

At Mach 10 (fig. 49), the range is of the order of 1100-1400 K. 

At Mach 15 (fig.50), small region situated near the leading edge have a temper- 
ature higher than 1800 K but none of them actually reached 1900 K. For this case, 
the temperature on the surface ranges from 1400 to 1700 K. 

At Mach 20 (fig.51), the situation is somewhat different. A thin region of laminar 
flow exists near the leading edge and the temperature drops in this region. But 
transition takes rapidly place after that region and the temperature rises above 
1800 K. 

The thermal protection proposed in chapter 3 section 6 seems to be sufficient to 
protect the lower waverider surface. Moreover, even in the most severe conditions, 
active cooling was not found to be necessary if a carbon-carbon protection can be 
expanded somewhat downstream of the leading edge. 

Finally, a last set of figures (52 to 55) presents the convective heat transfer and the 
wall temperature distribution along the trailing edge of the upper and lower surface of 
the vehicle. The main purpose here is to illustrate the difference in wall temperature 
on the upper and lower surfaces as well as the large temperature gradient from one 
streamline to the other. These large temperature gradients will result in heat transfer 
by conduction through the structure that will have to be managed by an appropriate 
disposition of heat sinks in the structure. 

Waverider performance 

The conventional viscous-optimized waveriders with a sharp leading edge, pro- 
duce high IVD values. Question : Does the finite leading edge bluntness gready 
affect the lift-to-drag ratio ? 
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Looking at Table 6, we see that the leading edge pressure drag remains a small 
fraction of the inviscid drag. Note that for the 8 cases presented, the leading edge 
drag represents 4 to 8% of the inviscid drag. Therefore, we can conclude that for the 
given carving law defined in chapter 2 section 3, the leading edge drag will represent 
a small but non-negligible fraction of the inviscid drag of the waverider. Note that 
Navier-Stokes calculations performed on blunted waveriders 6 have shown similar 
results in terms of the leading edge drag. Since the waverider considered presents 
similar geometry at a given Mach number, very small differences are observed in 
the inviscid drag and in the inviscid value of L/D. On the other hand, due to their 
differences in Reynolds number the viscous drag at a given Mach number is more 
affected. Although the high dynamic pressure is five times larger than the lower 
dynamic pressure, the unit Reynolds number based on the freestream conditions is 
more than five time larger at high dynamic pressure. Since the skin friction is inversely 
proportional to some function of the local Reynolds number, the viscous drag will 
be smaller at high dynamic pressure. At low dynamic pressure, the viscous drag 
decreases with increasing Mach number due to a progressive laminarisation of the 
vehicle surface. At high dynamic pressure, the viscous drag seems to be constant 
with Mach number. In fact, at Mach 5 and 10, the leading edge is fully turbulent 
and in both cases the waverider is wetted by turbulent flow. At Mach 15 and 20, the 
leading edge is laminar but the local Reynolds number is still so high that most of 
the waverider surface is wetted by turbulent flow. 

Comparing the pressure drag and the viscous drag, we see that the pressure drag 
is always higher than the viscous drag. Bowcutt^ and Corda 3 demonstrated that the 
highest L/D value is always obtained when the pressure drag and the viscous drag 
are approximately equal. This suggests that higher L/D values could be obtained 
for smaller cone angle in the 8 cases treated here. Nevertheless, when compared 
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to the best optimum in L/D obtained by Bowcutt 2 and Corda 3 (fig.56) we see that 
the present values are generally lower except for Mach 5, high dynamic pressure 
where the L/D obtained here is slightly higher than the new L/D “barrier” defined by 
Bowcutt 2 . At low dynamic pressure, where the entire leading edge was laminar from 
Mach 5 to Mach 20, the variation of L/D with Mach number (Table 6) presents a local 
minimum at Mach 15 similar to the results of Bowcutt 2 and Corda 3 (fig.56). At high 
dynamic pressure, this local minimum is not observed but the Reynolds numbers are 
not comparable. 

Waverider optimized for volume 

All the previous test cases discussed up to here have been optimized in order to 
maximize L/D. They all present the same type of triangular planform shape with an 
almost constant sweep angle from the second to the last point of the leading edge. 
If another objective function is chosen for the optimization algorithm, for instance 
waverider volume, the code comes out with a totally different planform shape. The 
test case at Mach 20 and high dynamic pressure was run with volume as the objective 
function. The leading edge results and the planform shape are given in figure 57 and 
58. As we see, a major difference appears in the sweep angle distribution. The 
sweep angle is no longer nearly constant but gradually increases along the leading 
edge. This has a major impact on the variation of R t . As a consequence, transition 
is now predicted near the nose and relam inarization is occurring later when R 0 drops 
below 250 (see chapter 3, section 2, topicl). Since turbulent flow is present along 
a large part of the leading edge, a relatively large portion of the upper surface will 
now be wetted by turbulent flow as opposed to the waverider optimized for L/D 
for the same flight condition which has a fully laminar upper surface. Even more 
surprising, for the same kind of planform shape, at low dynamic pressure the flow will 
be turbulent along the leading edge for the Mach 5 and Mach 10 cases. Therefore, 
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due simply to a change of planform shape, our previous observation about the nature 
of the flow for the LTD optimized waverider is totally reversed. In the case of the 
waveriders optimized for volume, only two test cases (Mach 15 and 20, low dynamic 
pressure) would present a fully laminar leading edge. This major difference will have 
an negative impact on the performance of the volume optimized waveriders. When 
compared to the maximum L/D (Table 7), the volume optimized wavenders , show 
a moderate decrease in L/D ranging from 8 to 15%. On the other hand, they show 
an increase in volume of the order of 60 to 350%. From the heat transfer point of 
view, the volume optimized cases present a higher average heat transfer rate due to 
smaller sweep angle and the turbulent nature of the flow (compare fig.35 and 37). 

Finally, due to the continuous increase of sweep angle along the leading edge, 
a favorable pressure gradient exists along the leading edge. As a consequence, 
the assumption of an infinite swept cylinder with no pressure gradient along the 
attachment line might not be even true locally anymore. What is the impact of a 
favorable pressure gradient on the stability of a 3D boundary layer at hypersonic 
speed ? The problem was partially studied by Zurigat and al. 22 . They concluded 
that the effect of pressure gradients on stabilizing three-dimensional first-mode waves 
is much larger than its effect on stabilizing two-dimensional first mode waves. But 
we know that at hypersonic speed, the second-mode waves are dominant. They also 
demonstrated that the effectiveness of the pressure gradient on the natural laminar 
flow control of 2D second-mode waves decreases at hypersonic speed but nothing 
is said about the effects of pressure gradients on a 3D second-mode waves. This 
problem is apparently not well understood 23 and our question will remain without 
answer for the present time. By lack of any useful information on favorable pressure 
gradient and a better model for our swept leading edge we will continue to locally 
assume an infinitely swept cylinder with no pressure gradient. 
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Chapter 5 Conclusions and recommendations 


The results of the present study lead to a number of important conclusions 
regarding the transition process along the attachment line of the leading edge and 
the heat transfer characteristc of a hypersonic waverider. 

Section 1 Conclusions 

A summary of the major results discussed in the previous chapter are itemized 
below : 

• Geometry : 

• It is possible on a 60 meter waverider to successfully carve a radius distri- 
bution along the leading edge without strongly affecting the planform shape. 
The radii obtained are sufficient to implement an active cooling device. 

• Transition at the attachment line of the leading edge : 

• for the L/D optimized waveriders, the flow was found to be laminar in all 
cases except at A/qo < 10 at high dynamic pressure. 

• for the volume optimized waveriders, the flow was found to be turbulent in 
all cases except at > 15 at low dynamic pressure. 

• Cooling requirements and thermal protection system : 

• for a maximum wall temperature of 1900 K and a minimum radius of 1 mm, 
active cooling is not required at Mach 5 for the altitudes considered. 

• for a maximum wall temperature of 1900 K and a minimum radius of 10 
mm, active cooling will be required for > 10 

• In the worst conditions (Mach 20, high dynamic pressure), the active cooling 
device proposed for the leading edge was always working within its cooling 
capabilities. 
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• In the worst conditions (Mach 20, high dynamic pressure) and for a cone angle 
of 7.5 degree , a ceramic composites thermal protection system is sufficient 
for the lower surface and no active cooling should be necessary. 

For a larger cone angle this might not be true. 

• Performance : 

• The inviscid drag due to the leading edge was found to be a small but non- 
negligible (4 to 8%) fraction of the inviscid drag. 

• Although the viscous drag at the leading edge can be neglected, the presence 
of the leading edge will greatly influence the transition front of the upper and 
lower surface and therefore affect the viscous drag of the entire vehicle. 

• When optimized for L/D, the performance obtained in this study are within 
the two L/D barriers defined by Bowcutt for similar flight conditions. 

• When optimized for volume, the geometry and the transition pattern are 
radically different than the waverider optimized for L/D. The decrease in 
L/D ranges from 8 to 15 % for an increase of 60 to 350 % in volume. 

Section 2 Recommendations for further study 

This study has always assumed that the flow was in a calorically perfect state 
at the wall. For a maximum wall temperature of 1900 K and a fully catalytic wall, 
this assumption is perfectly valid in a region of high pressure. In the low pressure 
regions this assumption is questionable since we know that low pressure will enhance 
dissociation. Therefore, it might be interesting to consider a study where pressure 
and temperature at the wall would be used to determine the amount of dissociated 
oxygen present on the vehicle surface. 

In another category, if we look to the geometry of the waverider optimized for 
L/D or maximized for volume both have to be considered as unrealistic. The first 
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ones for an evident lack of volume and the second ones for a bad distribution of 
this volume. The optimization algorithm creates high volume by keeping the base 
thickness almost constant and by increasing the span as much as it can. This results 
in a very flat waverider with large span. These configurations are probably useless 
because the structure needed to support this large span would be very heavy. An 
easy way to produce useful volume is to calculate the moment of inertia of a cross 
section from the symmetry plane. The moment of inertia of a cross section could 
be calculated at different stations along the waverider and averaged. The averaged 
moment of inertia has the dimension of m 3 and can be used to non dimensionalize 
the volume. If the following objective function is chosen : 



we should obtain vehicles with good values of L/D and a optimum distribution of 
this volume. L/D or the volume distribution could eventually be stressed according 
to the value of n. 

Finally, the present study does not show any results for a 10 meter waverider. 
For a vehicle of that size and the carving law defined by equation (1), it is not 
possible to implement everywhere along the leading edge the minimum radius of 10 
mm required by the active cooling device. By increasing the generating cone angle, a 
larger deflection of the streamlines can be obtained as well as larger values of the radii. 
But, due to the increased deflection of the streamlines, the performance of these kind 
of vehicles drops considerably and the wall temperature on the lower surface is much 
higher than the maximum temperature allowed by the thermal protection. Therefore, 
vehicles of that size will require another type of active cooling at the leading edge 
allowing smaller values of the leading edge radii and higher cooling performance. 
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Appendix A Tables 


Test cases 

M 0 o 

Qo o (atm) 

Altitude 

Reyoo * 10 6 

Vooimjs) 




(km) 

(1/m) 


1 

5 

0.196 

30.4 

1.625 

1529 

2 

10 

0.199 

40.4 

0.738 

3243 

3 

15 

0.203 

46.8 

0.458 

5047 

4 

20 

0.206 

51.5 

0.347 

6744 

5 

5 

0.998 

19.8 

9.383 

1475 

6 

10 

0.990 

28.8 

4.375 

3024 

7 

15 

0.998 

34.4 

2.691 

4699 

8 

20 

1.007 

38.6 

1.914 

6423 


Table 1 Flight conditions 


Test 

cases 

Moo 

Q oo 

(atm) 

Cone 

< 

L(m) 

1 

Rough 

Cooling 

1 

Max 

Box 

Min 

Sir 

1 

5 

0.2 

7.5 

60 

m 

Yes 

No 

1900 

0.4 

.075 

2 

10 

0.2 

7.5 

60 

1.4 

Yes 

Yes 

1900 

0.4 

.075 

3 

15 

0.2 

7.5 

60 

1.4 

Yes 

Yes 

1900 

0.4 

.075 

4 

20 

0.2 

7.5 

60 

1.4 

Yes 

Yes 

1900 

0.4 

.075 

5 

5 

1.0 

7.5 

60 

1.4 

Yes 

No 

1900 

0.4 

.075 

6 

10 

1.0 

7.5 

60 

1.4 

Yes 

Yes 

1900' 

0.4 

.075 

n 

15 

1.0 

7.5 

60 

1.4 

Yes 

Yes 

1900 

0.4 

.075 


20 

1.0 

7.5 

60 

1.4 

Yes 

Yes 

1900 

0.4 

.075 


Table 2 Input data files options 
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Test 

cases 

Surface 

emissivity 

Minimum 

Radius 

Rcycrit 

C(0) 

C(l) 

m 

1 

0.8 

1 mm 

15.10 6 

2% 

50% 

0.9 

2 

0.8 

1 0 mm 

Sheetz 

2% 

50% 

0.9 

3 

0.8 

1 0 mm 

Sheetz 

2% 

50% 

0.9 

4 

0.8 

1 0 mm 

Sheetz 

2% 

50% 

0.9 

5 

0.8 

1 mm 

15. 10 6 

2% 

50% 

0.9 

6 

0.8 

1 0 mm 

Sheetz 

2% 

50% 

0.9 

7 

0.8 

1 0 mm 

Sheetz 

2 % 

50% 

0.9 

8 

0.8 

1 0 mm 

Sheetz 

2 % 

50% 

0.9 


Tabic 3 Code options 


Test cases 

Shock 

Base/L 

Box 


Veff 


< 



■Bl 


1 

13.7 

0.077 

0.34 

1669 

0.105 

2 

9.9 

0.106 

0.22 

1845 

0.170 

3 

9.0 

0.116 

0.19 i 

1685 

0.193 

i 4 

8.7 

0.119 

0.19 

1850 

0.200 

5 

13.7 

0.075 

0.35 

1685 | 

0.103 

6 

9.9 

0.105 

0.23 

1931 

0.169 

7 

9.0 

0.117 

0.15 

1389 

0.207 

8 

8.7 

0.120 

0.13 

1093 

0.228 


Tabic 4 Waverider dimensions* LTD optimized 
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R. 

max 

V 

max 

(mm) 

9w 

max 

(kW/m 2 ) 

T 

1 w 

max 

(*) 

(%) of 
turbulent 
flow 

(%) of 
cooling 

272 

0.159 

73.1 ; 

1127.1 

0.0 

0.0 

223 

0.451 

842.9 

1900.0 

0.0 

0.1 

154 

0.701 

2314.1 

1900.0 

0.0 

3.3 

118 

0.873 

5067.4 

1900.0 

0.0 

100.0 

623 

0.067 

84.1 

1167.3 

92.0 

0.0 

538 

0.185 

1500.2 

1900.0 

91.0 

0.6 

304 

0.290 

5233.8 

1900.0 

0.0 

100.0 

250 

0.436 

9545.7 

1900.0 

0.0 

100.0 


Table 5 Leading edge results, L/D optimized 


Cdie 

Cd p 

L/D 

Cdyisc 

Cdtot. 

L/D 

*10 3 

*10 3 

invisc 

*10 3 

*10 3 


0.19 

2.6 

13.5 

2.1 

4.7 

mm 

0.24 

4.0 

8.7 

2.0 

6.0 

5.8 

0.27 


7.8 

1.6 

6.1 


0.23 

4.6 

WKEU 

1.4 

6.0 

5.8 

0.20 

2.5 

13.7 

1.4 

3.9 

8.7 

0.26 

4.0 

8.7 

1.4 

5.4 

d 

0.19 

4.5 

mm 

1.4 

5.9 

6.0 

0.25 


mam 

1.4 

5.9 

5.9 


Table 6 Wavcridcr performances, L/D optimized 
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Tests 

cases 

L/D 

optimised 

Volume 

(™ 3 ) 

LTD 

Volume 

optimised 

(>» 3 ) 

AWD 

(%) 

AVol 

(%) 

5 

8.7 

1685 


2763 

-8.0 

64.0 

6 

6.4 

1932 

nm 

4341 

-10.9 

124.8 

7 

6.0 

1388 


5271 

-15.0 

279.8 j 

8 

5.9 

1093 

5.0 

5089 

-15.3 

356.6 


Tabic 7 Comparison of waverider optimized for L/D and optimized for Volume 
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Appendix B Figures 




Figure 2 Illustration of five original leading edge points for a typical case (from reference 4) 
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Figure 3 Leading edge carving to accommodate bluntness 



Figure 4 Comparison of the original and the modified planform after carving 
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Figure 9 Normal cross section to the leading edge used for pressure integration 



Figure 10 The infinite swept cylinder 
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Figure 11 Flow near ihc leading edge of a swept cylinder (from reference 11) 



Figure 12 The variation of with d/ rf and s/r} for 
the appearance of first bursts of turbulence (from reference 12) 
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Figure 15 Co-ordinate system to investigate the instability 
of a three dimensional boundary layer (from reference 11) 



Figure 16 Variation of the velocity profile with the disturbance propagation angle (from reference 11) 
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Figure 17 Variation of the cross flow Reynolds number at 
transition onset with stream wise shape factor (from reference 11) 



SWEEP ANGLE (A) t d*0 


Figure 18 Effect of sweep angle on the normalized transition Reynolds number (from reference 2) 
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LEAST-SQUARES CURVE FIT 

Lo 9 10 (Re x ) t • 6.421 Exp (0.0001209 
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Figure 19 Correlation of transition Reynolds number with Mach number (from reference 2) 

• — o Bo wcutt's curve fit 
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q w v'l'n/Pi* 



Figure 21 Stagnation point heal transfer coefficient versus Knudsen 
number = 7.5 km/s, R„ = 1 inch (from reference 17) 



Figure 22 Stagnation point convective healing comparisons (from reference 17) 
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Carbon-carbon 
structure 



Spanwise heat pipes 
(stagnation region only) 
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H 2 coolant 


Heat pipes 
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Figure 23 Carbon-carbon / refractory metal heat-pipe concept for a wing leading edge (from reference 20) 
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Figure 24 Performance of leading edge cooling device (from reference 20) 
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Figure 25 Constant dynamic pressure trajectories 
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Mach = 5 


Mach r 10 



Mach = 15 


Mach s 20 



Figure 26 Waveriders optimized for L/D, Qoo = 0.2 atm. 



Figure 27 Waveriders optimi/cd fm L/D, Q ^ = 10 atm 
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Figure 31 Leading edge results. L/D optimized. A/., - 10. - 10 aim 
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Figure 32 Leading edge results, L/D optimized, M «, = 15, Q ^ = 0.2 atm . 
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Figure 33 Leading edge results, L/D optimized, M <>. = lf>, Qc*. — 10 atm. 
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Figure 35 Leading edge results, L/D optimized, — 20, 


1.0 aim. 
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Flow nature and grid details 



Figure 36 Transition front, M <*, — 5, Q v , = 0.2 atm. 
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Flow nature and grid details 



Figure 37 Transition front, M = 5 , Q& = 1-0 atm. 




Flow nature and grid details 
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Figure 38 Transition front, A/^ = 10, Q = 0.2 atm. 
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Figure 39 Transition front, M x , = 10, Q = 1.0 atm. 
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Flow nature and grid details 

Planform view 
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Flow nature and grid details 



Figure 41 Traasilion from, M^, — 15, Q. x = 10 “l" 1 
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Flow nature and grid details 
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Figure 42 Transition front, = 20, Q ^ = 0.2 atm 
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Flow nature and grid details 



Figure 43 Transition front, Afo,. = <J0, Qot — 1-0 aim 
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Convective heat transfer on the upper surface 
along different streamlines 



Convective heat transfer on the lower surface 



1.00 X/L 


Figure 44 Convective heat transfer on the upper and lower surface 
along different streamlines, = 5, = 1.0 otm. 
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Convective heat transfer on the upper surface 
along different streamlines 


(kW/M 1 ) Machs 10.0, 79 t. 0.991, RaytC* •)« 



Convective heat transfer on the lower surface 
along different streamlines 



X/L 


X/L 


Figure 45 Convective heat transfer on the upper and lower surface 
along different streamlines, Moo = Q <*> = 1 0 afm. 
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Convective heat transfer on the upper surface 
along different streamlines 

(kW/M*) Mach. 15.0, AK(km). 34.4, 0<atm). 0.001. RayL(‘a-<)« 111.5 



X/L 


Convective heat transfer on the lower surface 
along different streamlines 

(kW/M 2 ) Mach. 15.0 . Ah (km)* 34,4, Q(alm)» 0.000, Rayl{*a-I)s 111. 5 



Rough leading adge, with krtamal cooling 


Figure 46 Convective heat transfer on the upper and lower surface 
along different streamlines, M = 15, Q <*, ~ 1.0 aim 
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Convective heat transfer on the upper surface 
along different streamlines 

Macfe 20.0, Ah(hirg« )l.l. 1.007, RayL(*a-0)« 114.0 


(kW/M 2 ) 



Rough Wading adg«, with Imamal cooling 


Convective heat transfer on the lower surface 
along different streamlines 

(kW/M 2 ) Machs 20.0, All(hm)« tt.l, Q(atm)« 1.007, Rayl<*a-0)« 114.0 



ooo 0.25 0.50 0.75 1-00 X/L 

Rough Wading adga, with Internal cooling 

Figure 47 Convective heat transfer on the upper and lower surface 
along different streamlines. Moo — 20, Q<^ = 1.0 atm. 
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Lower Surface 



Figure 48 Wall temperature distribution on the lower surface. Moo = 5, Qqc = 1.0 atm 
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Lower Surface 
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Lower Surface 



i- 


Figure 51 Wall temperature distribution on the lower surface, = 20, Qoo — 1-0 atm. 
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Convective heat transfer along the trailing edge 


(kW/M 3 ) Mach. S O, Ak(km). 1».i, Q(atm}« O.fOS, RayLTal). 9S3.0 



canlarltna 


Rough laading tdga, no Internal coding 


wing lip 


Y/L 
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Figure 52 Convective heal transfer and wall temperature along the trailing 
edge of the upper and lower surface, Moo = 5, Qoo = 1.0 atm. 
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Figure 56 Wavcriders performance comparison 
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Figure 57 Leading edge results, Volume optimized, Mqc = 20, Qoo = 1.0 atm. 
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Flow nature and grid details 



Figure 58 Transition front. Volume optimized, M ^ = 20, Qoo = 1-0 aim. 
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Figure 59 Waverider optimized for L/D Vs Waveridei optimized for Volume, 
Moo = 20, Q«, = 1.0 atm. 
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Appendix C Constraints on waverider 
construction and simplex optimization 


The following section was almost entirely copied from McLaughlin 4 . We have 
simply added here the set of new constraints which are now imposed in the modified 
version of the code. 

Certain constraints must be enforced in order to create valid waveriders, and to 
have control over the optimization procedure. Each constraint is stated and discussed 
in the following sections. 

Recall that the leading edge curve projected in the x-y (fig.2) plane is the initial 
parameter that needs to be specified in order to create a single waverider. Past 
experience in generating waveriders, gained by Bowcutt and Corda, has shown that 
there are a number of restrictions that need to be placed on the geometry of each 
leading edge in order to insure that a valid waverider is constructed. The reasons for 
some of these constraints are relatively obvious, while others are somewhat subtle. 

Original constraints set by Bowcutt, Corda, McLaughlin : 

Constraints on the leading edge shapes: 

1. The second leading edge point away from the centerline cannot be on the line of 
symmetry. — If this were to occur, the resulting waverider would have a sharp 
spike on the upper surface which is undesirable. In addition, two streamlines 
would lie in the same flow plane. This, too should be avoided because problems 
will arise in the integration of the waverider volume and surface areas. 

2. All leading edge points must lie within the shock wave. — Since the flowfield is 
not defined outside of the shock wave, the streamlines emanating from these points 
would have to be traced as freestream streamlines, resulting in extra surfaces. 
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More importantly, the shock wave would not be attached to the leading edge, 

which violates the basic waverider premise. 

3. The radial distance to each of the leading edge points from the centerline to 
the shock wave must continually increase. - In some cases, a violation of this 
occurrence will cause the waverider’ s upper and lower surfaces to cross over each 

other. This is obviously an unrealistic result. 

4. The polar angle to each of the leading edge points from the centerline to the 
shock wave must continually increase. — A violation of this occurrence also 
causes the waverider’ s upper and lower surfaces to cross each other. 

A completely different set of constraints are placed on the entire vehicle geometry 
for the optimization procedure. For each of these cases, even though a valid leading 
edge is created, the resulting waverider may exceed other limits imposed by the user. 
The user imposed limits provide a control mechanism over the optimization procedure. 
The simplex is never allowed to move beyond these bounds, and therefore some 
control is gained over which local minimum is found by the optimization procedure. 
All of the constraints listed below have the option of being made inactive. This is 
sometimes done since it is occasionally necessary to relax some of the constraints to 
allow the optimization procedure to continue operating properly. 

Constraints on the Vehicle Geometry: 

1 . Minimum Slenderness Ratio - The user can adjust this constraint to limit the 
lower bound of the slenderness ratio of the vehicle. Any vehicle that has a 
slenderness below this minimum is discarded by the optimization procedure. 

2. Box Size — The user can set both the upper and lower bounds on the box size of 
the vehicle. The box size is defined as the semi-span to length ratio of the vehicle. 
This constraint can be imposed to keep the waverider planform narrow, or it can 
be relaxed to allow the waverider to have a wide planform. The suggested values 
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for these parameters are BOXMAX=0.4 and BOXMIN=0.1. These are the values 
that arc used for most of the applications in the present study. Any vehicles 
outside of these bounds arc discarded by the optimization procedure. 

3. Volume — The user can also impose limits on the minimum and maximum 
volume of the waveriders. The minimum limit allows the user to insure that none 
of the waveriders fall below a given limit. This is useful when optimizing for 
maximum L/D since no waverider will be considered if it has a volume below 
this minimum. The maximum limit is useful when relaxing either the slenderness 
ratio or box size constraints. This keeps the waveriders from getting drastically 
out of proportion in some of these cases. 

New set of constraints defined in the present study : 

Constraints on the Vehicle: 

1. Minimum radius : When cooling is needed, the value of the minimum radius 
is set to 10 mm. When cooling is not needed the minimum radius is set to 1 
mm. Note that in both cases, the minimum radius needed to avoid low density 
effects on the heat transfer is also calculated. Therefore, the code will always 
consider the maximum value of these minimum radius as the final constraint on 
the radius distribution. 

2. Planform shape : In order to avoid unrealistic shapes, the waverider will be 
rejected if the sweep angle distribution presents a local minimum. 

3. Maximum wall temperature : If the maximum temperature is exceeded, when 
cooling is not pre-required, the waverider will be rejected. If cooling is pre- 
required and the maximum wall temperature is exceeded, the program will check 
that the cooling requirement is within the possibility of the active cooling device. 
If not, the waverider is rejected. 
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study shows that surface radiative cooling might be sufficient. In all cases, 
radiative cooling is sufficient for the upper and low surfaces of the vehicle if 
ceramic composites is used as thermal protection. 
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